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ABSTRACT 



Aims. We uncover the Ht flows in the Corona Australis molecular cloud and in particular identify the flows from the Coronet cluster 
Methods. A deep, near-infrared H2 v=l-0 S(l), 2.122/jm-line, narrow-band imaging survey of the R CrA cloud core was carried 
out. The nature of all candidate-driving sources in the region was evaluated using data available from the literature and also by fitting 
the spectral energy distributions (SED) of each source either with an extincted photosphere or YSO model. Archival Spitzer-IRAC 
and MIPS data was used to obtain photometry, which was combined with USNO, 2MASS catalogs and millimeter photometry from 
the literature, to build the SEDs. We identify the best candidate-driving source for each outflow by comparing the flow properties, 
available proper motions, and the known/estimated properties of the driving sources. We also adopted the thumbrule of outflow power 
as proportional to source luminosity and inversely proportional to the source age to reach a consensus. 

Results. Continuum-subtracted, narrow-band images reveal several new Molecular Hydrogen emission-line Objects (MHOs). 
Together with previously known MHOs and Herbig-Haro objects we catalog at least 14 individual flow components of which 11 
appear to be driven by the R CrA aggregate members. The flows originating in the Coronet cluster have lengths of ~0. 1-0.2 pc. Eight 
out of nine submillimeter cores mapped in the Coronet cluster region display embedded stars driving an outflow component. Roughly 
80% of the youngest objects in the Coronet are associated with outflows. The MHO flows to the west of the Coronet display lobes 
moving to the west and vice-versa, resulting in nondetections of the counter lobe in our deep imaging. We speculate that these coun- 
terflows may be experiencing a stunting effect in penetrating the dense central core. 

Conclusions. Although this work has reduced the ambiguities for many flows in the Coronet region, one of the brightest H2 fea- 
ture (MHO2014) and a few fainter features in the region remain unassociated with a clear driving source. The flows from Coronet, 
therefore, continue to be interesting targets for future studies. 

Key words. Stars:formation 



1. Introduction 

The Corona Australis molecular cloud (see the review by 
iNeuhauser & Forbich 2008) . located in the southern hemisphere 
of the sky, is one of the nearest (130-170pc; recommended value 
is 130pc) and best studied star-forming regions with intense 
star formati on activity. The 1. 2mm dust continuum emission, 
mapped by IChini et al. I (l2003h using the SIMBA array on the 
SEST telescope, provides a good representation of the dense 
regions of this cloud where most of the active star formation 
takes place. The densest region of this cloud roughly coincides 
with the Herbig Ae/Be star R Corona Austrinae (RCrA), sur- 
rounded bj;_acluster_ofyoui^^ well-known as the Coronet 
cluster dTavlor & Storev 1 11984|) . The Coronet cluster and as- 
sociated molecular core has been the target of several studies 
at different wavelengths, from X-rays (for a complete list, see 
the review bv N euhauser & Forb ich 2008) through the optical 
and infrared (IWilking et al. to the submillimeter (submm) 

(iNutter et al. Il2005h . 

Outflows and jets are ubiquitous phenomena associated with 
star formation activity. Bipolar outflows are unambiguously 
identified by mapping the two lobes using molecular line transi- 
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tions from species such as CO or SiO at (sub)mm wavelengths. 
Observations of optical Herbig-Haro (HH) objects and/or near- 
infrared (near-IR) H2 v=l-0 S(l) Une emission features often 
complement these radio data; indeed, they can be more success- 
ful in tracing the flow components, particularly in crowded re- 
gions, because of the higher spatial resolution typically available 
at these shorter wavelengths. 

Molecular line maps of the Coronet region show a promi- 
nent east -west bip olar outflow and a weak north-south flow 
(Andersoi ietal. 11 997: Groppi et al. 2004, 2007). The east-west 
flow is thought to be driven by sources around IRS7, which is 
located roughly 15-20" southeast of RCrA. The weak north- 
south flow is driven by the IRSl source. This outflow is asso- 
ciated with HH 100, and its driving source is often referred to 
as HHIOO-IR. In contrast, at least 20 HH objects are known 
to be located in t he close vicinity of the Coronet cluster (see 
IWang et ani2004 and references therein). These data suggest 
the presence of several flow components. The east-west molecu- 
lar outflow has been mapped in a number of molecular lines (e.g. 
CO, HCO-t-, and SiO), and these observations li kewise infer the 
prese n ce of more than one outflow component (lAnderson et al. I 
119971: iGroDDi et al. I I2OO4I) . However, decomposition of these 
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flow components has not been possible using the lower spatial 
resolution radio data. 

The optical observations of IWang et al. I ( 2004 ) and, more 
recent ly, th e near- infrared H2 ima ges of ICaratti o Garatti et al. I 
(1200 6^ and ' Peterson et al. I (1201 ll) provide a far superior view 
of outflow activity in the region. Even so, the concentra- 
tion of candidate driving sources within a small region in the 
Coronet cluster, together with the large number of flow com- 
ponents, has made it difficult to unambiguously associate the 
flows with their driving sources. The rather complete census of 
the embedded populati on of young stars in the re gion uncov- 
ered by the new X-ray (Forbrich & Preibish '2007^ and Spitzer 
Space Telescope obs ervations (C unle_&.Sidlia-.A g.uilai:., 201 1; 
iPeterson et al. II20TTI) has prompted a fresh look at this region. 

In an attempt to effectively trace the outflow activity 
in this region, we have conducted a deep H2 v=l-0 S(l) 
2.12jumnarrow-band imaging survey of the northern part of 
the Corona Australis cloud. Many new Molec ular Hydrogen 
emission-line Objects (MHOs lOavis et ani2010h are discovered 
from the new data, increasing the number of known flow com- 
ponents. In this work, we present and organize all the new and 
previously known outflow components. 

For this purpose we use the following strategy. It has been 
suggested that the outflow power (and therefore MHO activity) 
is directly proportional to the luminosity and inversel y to the 
evolutionary stage of their driving sources (e.g Bontemp s et al, I 
119961) . First, by using the available information in the literature 
and then by modeling the spectral energy distributions (SED) of 
candidate driving sources in the region, we obtain a handle on the 
luminosity (mass) and evolutionary stage (age). Combining this 
information with the known properties of the outflows, as stated 
above, we then disentangle the individual flows and attribute the 
best driving sources to each flow. We also use the information 
on the proper m otions of many knots in this region presented by 
iPeterson et al. I (1201 ih to arrive at the best-case scenario. 



2. Observations and data analysis 

2. 1 . Narrow-band H2 imaging 

Near-infrared observations of the Corona Australis molecular 
cloud were carried out using the 4m Vitor Blanco Telescope 
at the Cerro-Tololo Inter-American Observatory (CTIO) on the 
nights of 11 to 13 July 2009. Imaging observations in the H2 
2. 12/zm narrow-band filter and in the 2. 2//m broad-band K fil- 
ter were obtained. These observations cover the northern part 
of dense g a s regi ons traced by the 1.2mm observations of 
IChini et al. I (l2003h . The near-infrared wide-field camera ISPI, 
used to obtain the observational data, has a 10' field of view and 
a plate scale of 0.3" per pixel. A nine-point jitter pattern was 
used to obtain the observations, and the jitter sequence was re- 
peated three times at each pointing of the telescope in each filter 
An exposure of 60 sec per jitter position was used for the H2 
observations, resulting in a total integration time of 27 mins. For 
the K-band observations, the exposure time per jitter position 
was kept at 5 sec, resulting in an integration time of 2 min. 

Each set of jittered images were median combined to obtain 
a skyframe that was subtracted from each individual image. The 
sky-subtracted images were subsequently flattened by dividing 
by a smooth surface-fit obtained from the image. The individual 
exposures were co-added and mosaicked together using KAPPA 
and CCDPACK routines in the StarlinlQ suite of software. 
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2.2. Spitzer pliotometry 

Archival Spitzer imaging data were used to obtain photometry 
of the sources in the Coronet region, which in turn were used 
for SED modeling. The data used correspond to the astronomi- 
cal observational requests (AOR's) 3650816 and 17672960 (PI: 
G. Fazio) and 27041280 (PI: L. Allen). The BCD data in each 
of the four IRAC channels were corrected using artifact mitiga- 
tion software provided by the Spitzer Science Center (SSC), and 
were then mosaicked together to obtain images corresponding 
to the region covered by our H2 narrow-band observations. We 
used the APP//(9r package in IRAF to obtain aperture photome- 
try using an aperture of 2.4"and a sky annulus of width 12" . The 
aperture corrections used for the Chl-4 IRAC bands are 1.213, 
1.234, 1.379, and 1.584. 

MIPS lAjjm data of the same region were obtained from 
AOR 3664640 (PI: G. Fazio). The BCD data were mosaicked us- 
ing standard templates available on MOPEX software. Aperture 
photometry was performed by using the MOPEX/APEX single 
frame template. Aperture size of 6", sky annuli of width 7"and 
an aperture correction factor of 2.05 were used. 



3. Spectral energy distributions and modeling 

As mentioned in the last paragraph of Sec.l, we will first ob- 
tain the information on luminosities (masses) and evolutionary 
states (ages) of candidate driving sources in the Coronet region. 
For this purpose, we first searched the literature to obtain the 
information derived from standard methods. When such infor- 
mation was not available, particularly for the candidate sources 
that are well-detected in the Spitzer observations, we analyzed 
their spectral energy distribution by modeling. By SED fitting, 
we first eliminated sources that display redenned photospheres 
and then used the SED modeling to make a relative comparision 
of YSOs based on the fitted parameters (class, luminosity, and 
accretion rates) to select the most probable driving source for 
each identified outflow. 

To this end we first built the SEDs of sources that appear 
to be associated with outflows and mod el these with a popular 
SED fitting tool dRobitaiUe et al. Il2007h . Not afl SEDs are uni- 
formly populated with the same photometric data because many 
of the millimeter sources are too bright (saturated) in the infrared 
bands. Some sources are resolved by optical and/or near-IR data 
and unresolved in the millimeter bands. Also, while we mod- 
eled all candidate driving sources in the mapped region, our 
focus is on a selected group of sources, which appear to drive 
outflows (see Table. 2). These sources are prominently infrared- 
millimeter visible and were all modeled with 4-6 "data points" 
and with 1-6 upper or lower limits. The photometric data, aper- 
tures used, and the source of the data are listed in Table. 1 . Only 
the IRAC and MIPS photometry are independently derived by 
us, and the remaining photometric fluxes and errors are directly 
obtained from the references listed in Table. 1. The apertures 
listed in Table 1 are those used for SED fitting and not nec- 
essarily the one with which photometry was obtained. This is 
because the SED fitter requires that the modeled source should 
not be larger than the aperture used for fitting. If a source is sat- 
urated in any band, the flux is assumed to be a lower limit and 
if the source is too faint or larger than the aperture used, then 
the flux is assumed to be an upper limit. The 2MASS magni- 
tudes are used as an upper limit if either the confusion flag was 
not "0" (implying confusion or extended source) or if the pho- 
tometric quality flag was not the best (anything other than A). 
Ten percent errors are assumed for fluxes for the shorter wave- 
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Table 1. Data used for SED fitting 



Band 


Aperture 


Reference 


USNO B, R, I 


3" 




2MASS J, H, K 


2.4" 




Spitzer IRAC 


3" 


this work 


SpitzerMWS lApm 


6" 


this work 


Spitzer MIPS 10pm 


9" 


this work 


450ytim SCUBA 


8" 


Nutter etal. (2005') 


850yum SCUBA 


14" 


Nutter etal. (2005) 


1.2mm SIMBA 


24" 


Chinietal. ('2003) 


1.1mm SMA 


3" 


GroDDietal. (2007) 



length bands (up to %pm ) and 20% errors were used for all other 
longer wavelength data. In most cases, the photometric accuracy 
is better than the assumed error. 

The resulting SED data were fed into the online SED fit- 
ting tool described by Robitaille et al. (2007). The distance and 
extinction values were allowed to vary in the range 130-170pc 
and 5^5 Al, respectively. The SED fitting tool attempts to fit 
the observed fluxes, either with an extincted photosphere or an 
young stellar object (YSO). Models are constrained by the good- 
ness of fit, x^- The YSO models used were obtained from a 
grid of 200,000 models, as described by Robitaille et al. (2006). 
Following the recommendations by .Robitaille et al. . (i2007i) . 
only models that satisfy the cnX.tn?iX^ best ~X^'^^ P^r data point, 
where x^best is for the best-fit model, are used in estimating rep- 
resentative values of some physical parameters. These are es- 
sentially the group of models displayed in Fig. 4. The results 
using this group of models are listed in Table. 3. The minimum, 
maximum, and weighted means (weights being the inverse of 
X^) are listed. It should be noted that the weighted means repre- 
sent the parameter value distribution better than the difference in 
the maximum and minimum values. This is because even a sin- 
gle outlier model or parameter within the chosen x^ criteria can 
produce a minimum and a maximum that are not representative 
of the overall distribution. The minimum and maximum values 
nevertheless represent the extreme values for a given parame- 
ter from the group of models shown in Fig. 4. The evolutionary 
stage of Class O/I quoted in Table 3 represents a fitted source age 
in the range l-10xlO"*yrs. 

In what follows we present the results from this procedure 
only for a few selected sources, for which the data available 
in the literature were not decisive in associating with detected 
MHOs. SEDs that were fitted with extincted photospheres are 
also not shown, except for one example, IRS 11 (see Fig. 4), as 
they were ruled out of qualifying as driving sources. 



4. The molecular yydrogen emission-line objects 
(MHOs) 

Figure 1 (online only) is a color composite generated by using 
the Spitzer IRAC-4.5//m image as red, the H2 narrow-band im- 
age as green and the SII narrow-band image (mosaicked using 
archival ESO data from Wang et al. (2004) as blue. The SII data 
do not cover the lower portion of the image. The scaling of the 
individual images is adjusted to effectively display the embedded 
sources and the outflow components simulatneously. Since [SII] 
and H2 trace different gas excitation conditions, they are often 
hugely complementary when one is trying to trace large sections 
of a protostellar outflow and/or connect outflow features with 
distant driving sources. 



The continuum subtracted H2 image of the mapped region 
is displayed in Figure 2 where five groups of MHOs are marked 
by ellipses. The majority of the MHOs are contained within the 
large central ellipse. In Table 2 we list the positions of the MHOs 
and the identified driving sources, along with the measured flow 
lengths in arcseconds. Of these, the newly detected MHOs 2005, 
2007, 2014, and 2017 are located away from the Coronet cluster 
and the dense core associated with it. The MHO features in the 
central region, corresponding to the Coronet cluster, are shown 
in Figure. 3. 

The Coronet region is composed of many previously 
known infrared sources plus new infrared coun t erparts 
for the millimeter sources ( Forbrich & Preibishi l2007t 



ICurrie & Sicilia-Aguilar]l201 it iPeterson et al. 11201 ih . Notably 
SMM4, SMAl, and SMM2 have been recently detected by 
Spitzer images. This demonstrates that these submm cores are 
protostellar in nature. Furtherm ore, the SMA observations by 
iGroppi et all (12007) and Peterson et al. I (12011) that provide 
millimeter fluxes are very useful in identifying the nature of 
the youngest protostars in the Coronet region. To disentangle 
each MHO flow and associate it with the best candidate driving 
source, we used the results of our SED modeling of all Spitzer- 
IRAC detected point sources (shown using circular symbols in 
Fig ure 3) and util i zed th e proper motion information provided 
by Peterson et al. I (1201 1 ). Together, these data have significantly 
improved the association of driving sources to the detected 
MHOs in the Coronet region. In Figure 4 we show the observed 
SED data and the fitted models for candidate driving sources. 
The results of our SED modeling are presented in Table 3. In 
the following, we first discuss MHO features that have driving 
sources reasonably well identified. 



4. 1 . MHOs with identified driving sources 

MHO 2000/HH 99 - SMA2: MHO2000 is the brightest, well- 
collimated bow shock feature to the north of Coronet, associ- 
ated with HH99. The suggestions of Wang et alj (|2004|) that it 
is driven by IRS6, and of D ayj.s et al. ( 1999 ) that it is driven 
by IRS9 were ruled out by Caratti o Garatti et aTl (|2006) based 
on the evolved nature of IRS6 and also on the po sition angle 
of proper motion vectors. Recentlv. iPeterson et al. I (2011) sug- 
gest that it is driven by IRS7 or CX034. The position angle of 
the proper motion vector f or MHO2000 was fou nd to be 39+6° 
bv ICaratti o Garatti et alTI ([2006) and 43+6° bv IPeterson et al. I 
(l2011h . These position angles of the proper motion vectors 
clearly point towards the sources SMA2, IRS7, and CX034. 

Since MHO2000 is a bright and well-collimated feature, we 
assume that it has to be driven by a ClassO/I type object in the 
region. SMA2 is bright at millimeter wavelengths, very bright 
(nearly saturated) at 24pm, and not detected at ^pm IRAC 
images, therefore, has insufficient data points to be modeled. 
Nevertheless, the above characteristics ar e strong enough to 
classify it as a Class protostar (see also iGroppi et al. 1 120071 : 
IPeterson et al. II20T1I) . In contrast, IRS7A display characteristics 
of being more evolved than SMA2 since it is bright at Sjum, 
not detected at (sub)millimeter, and also modeled to have a 
very low disk mass (0.004 Mq). Similarly, CX034 is modeled 
as a substellar object with very low disk mass as well. In 
fact, using the SMA n ondetections (both IRS7A and CX034), 
Peterson et al. I (1201 ih suggest an upper limit of only a few 
Jupiter masses for the disk in these sources. On the contrary, 
they estimate a disk mass of 0.063 Mq for SMA2, higher by 
a factor of ten. Therefore, we conclude that SMA2 drives the 
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MHO2000 flow, which is a bright, wen-collimated bowshock. iMever & WiUdngI (l2009h adopt a spectral type of FO for TCrA. 



MHO 2003,2008 (2007?) -IRSI/HHIOOIR: The source 
IRSl, also known as HH lOOIR is k nown to drive th e 
Herbig-Haro objects HH 100 and HH 97 ^miEK\^^- 
ICaratti o Garatti et all (l2006l) discovered MHO 2003 (an arc 
and an elongated knot) associated with IRS 1 . The deeper images 
presented here and by Pete rson et al. (201 li) reveal a new knot 
named as MHO2008, which is probably the tip of the flow 
associated with IRSl. While MHO2003 and MHO2008 are 
certainly driven by IRS 1 , it is likely that this flow is much longer 
that it terminates in MHO2007 (see the description of this MHO 
in the next subsection). MHO2003 appears to have a wide 
opening angle, MHO2009 from IRS2 and the confused flows 
from RCrA being the only two oth er cases in the region with 
large opening angles. Nisi ni et aTl (2005) estimate a spectral 
type of K-M for the star, however, they point out that IRS 1 and 
IRS2 are the only two sources with a bolometric luminosity 
that is 60-80% higher than the stellar luminosities. Also, using 
Brackett line characteristics, they suggest it probably originates 
in an expanding stellar wind rather than in accreting regions. 
The arc-shaped MHO feature 2003 may therefore be due to 
fluorescence rather than due to a large opening angle of the flow. 

MHO 2005 C, D, E-S CrA: This group of MHO knots 
represents the H2 counterpart of the well-kn own optical feature s 
HH 82A, HH 82B, and HH 729A,B,C dWang et all I200I . 
which are located near the bright T-Tauri star S CrA. This 
star coincides with the millimeter peak MMSl described by 
IChini et aT 1 120031) . Previous studies have suggested that these 
HH objects are dr iven by S CrA (.Rei purth & Graham 1988 ; 
IWang et alTI l2004. and so are these MHO knots. The knots 
MHO2005A B ('P eterson et al. ll20Tlb are not associated with S 
CrA, and is probably driven by IRS5a/b as discussed in Sect 4.2. 

MHO 2006/HH 7 30-IRS2: This source, first identified by 
IWilking et al. I (Il997h. is cla s sified a s a young star with a mass 



of I.4M0 by Nisini et al. ( 20051) . The recently discovered 



MHO 2006 feature (Peterso n et al. I I2OII) is clearly associated 
with this source and has a wide opening angle, most likely due 
to reasons like those discussed above with respect to IRSl. 
The HH objects HH 730A, HH 730B, and HH 730C, which 
terminate in a well-defined bowshock, clearly align with this 
source. We therefore associate HH 730 and the newly detected 
MHO 2006 object with a single flow driven by IRS2. 

MHO 2011 A, B & 2012 - SMM2: The newly discovered 
knots 2011 A, B, and the 2012 are located to the west of the 
Coronet cluster It comprises a number of knots and arcs that 
probably form part of at least two bow shocks. This flow is 
well aligned with the embedded IR source associated with 
the submm peak SMM2. The SED modeling presented here 
indicates that the Spitzer detected infrared counterpart to SMM2 
is a young low mass star of Class O/I type (see Table. 3). 

M HO 201 3/20 1 5/HH733 - TCrA: The faint bowshock MHO 
2013 dPavis et al . 1999) by the side of MHO20( )0 is the H2 
counterpart of the Herbig-Haro object HH 733. I Wang et alTI 
(12004 suggest that this flow is driven by TCrA. The newly 
discovered MHO 2015 is a clear bow-shock feature, lying 
to the south of TCrA, and it marks the southern lobe of the 
bipolar outflow originating from TCrA. This is probably the 
only unambiguously detected bipolar outflow traced by two 
complementing bow-shock features in the entire Coronet region. 



MHO2008D - Sta rl-100: The star identified as 1-lOOHH 
by IWang et aTl ( |2004|) is associated with a ne wly discovered 
faint H? knot that is labeled as MHO2008D by dPeterson et al. I 
I20T1I) . This is very likely an independent flow driven by this star. 

MHO 2016: This MHO feature has two conical nebulae A 
and B centered on TY CrA and HD176386, respectively. MHO 
2012A is an extended conical nebula, centered on the H^ emis- 
sion star TY CrA. MHO 2012B is a relatively faint, localized 
emission on the binary pre-main-sequence star HD176386. 



4.2. MHOs with ambiguous driving source associations 

MHO 2001 - IRS6b, CX034, IRS5N: MHO 2001 is one of the 
brightest features to the east of Coronet cluster, well immersed 
in the reflection nebula. There is much confusion in this area. 
Since it is one of the brightest MHO features, it should be driven 
by a relatively young object of Class O/I type. However, it ap- 
pears to be close to most of the can didate driving sources ( unlike 
MHO2000 which is farther awav). iPeterson et al. I (1201 lb quote 
CX034 as a possible candidate driving this flow. However, the 
same authors quote a position angle of the proper motion vector 
at 81° for the bright knot MHO2001A, which when traced 
backwards, points more directly to IRS6 or roughly towards 
IRS5N (considering an average of the PA from all the knots as- 
sociated with MHO2011). IRS6 is known to b e a binary sourc e 
with one young and one old component (iNisini et al. I l2005h . 
Following the proper-motion vectors, a likely scenario is that 
the younger object of this binary component drives MHO2001. 
Alternatively, MHO2001 could be driven by SMM4/IRS5N, 
which is a millimeter peak associated with a faint IR source, that 
is modeled as a very young substellar object. We also point out 
the low flux density of the SM A peak associated with IRS5N in 
comparison to SMAl or SMA2 lPeterson et al. I ( 1201 Ih . We favor 
the scenario of MHO2001 driven by IRS6 rather than IRS5N 
because of the likely stunting effect discussed in Sect.5 

MHO 2002, 2004 - RCrA : RCrA is the principal component 
of the Coronet cluster. We associate the previously detected 
bow-shaped object M HO 2002 with RCrA (labeled "E" by 
,Caratti o Gar atti et al. 1 12006) and possibly a knot from within 
the group labeled MHO 2004. Together, these would delineate 
an east-west bipolar outflow. These knots appear as a bowshock 
from a relatively wide-angled wind. The morphology of these 
MHO features suggests a wide-angled, poorly collimated wind 
rather than a collimated outflow. RCrA is a bright source that is 
saturated in all of the Spitzer bands. These properties conform 
to the nature of RCrA as a Herbig Ae/Be star with an estimated 
mass of 5-10Mo(Tuthin et al. 2001). Recently, Kraus et alT] 
(2009) have used the Very Large Array (VLA) Interferometer to 
study this source and find a disk that is oriented north-south, 
suggesting an east-west outflow. The proposed MHO associa- 
tions with RCrA presented above are in good agreement with 
this scenario. 

MHO 2004/HH 104 - IRS6b or RCrA: We associate HH 104 
and MHO 2004 with a faint, jet-like feature that is located close 
to the source IRS6. ICaratti o Garatti et al. I (l2006h suggest that 
IRS6 is too evolved to be driving a powerful outflow. In their 
high-resolution VLT data, iNisiniet al.1 (l2005h identify IRS6 as 
a binary system (IRS6a and IRS6b), and by using photometric 
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data, show that IRS6b has a stronger H-K color than IRS6a. Our 
K band data also partially resolve this binary source. Therefore, 
the HH 104/MHO 2004 outflow is probably driven b y the fainter 
- thou gh redder - infrared source, dubbed IRS6b by Nisini et al. I 
(l200l . IRS6 is unresolved in the mid-IR and submm photomet- 
ric data. 

MHO2004 also contains a bright bow-like knot that indi- 
cates a proper motion vector moving to the west; therefore, it 
could just as well be due to RCrA flow. Given that RCrA is 
a Herbig Ae/Be type star, the outflow from such a star would 
be composed of strong stellar winds and contribute to poorly 
coUimated flows composed of many knots and shock surfaces. 

MHO 2005 A, B - IRSSa/b: The bright infrared source 
IRS5 (iWilking et al. I [ 1997) is found to be a subarcsecond 
binary source dChen & Graham 1993; Nisini et al. 2005). 
The main component is labeled IRS5a by iNisini et al. I (l2005h 
and classified as a Class I source with low accretion activity. 
The faint, though well-defined, elongated MHO feature MHO 
2016 (Fig. 3) that terminates in a small bow shock is aligned 
well with IRS5a. The Class I na ture and the low a ccretion 
activity attributed to this source bv iNisini et al. 1 d2005h IS con- 
sistent with the low b rightness and morphology of MHO2016. 
iPeterson et al. I (1201 Ih associate this f eature wifli the pro minent 
east- west CO flow. Although IRS 11 (Wilking e ta"nil997 ) is lo- 
cated along the same projected line as MHO 2016, it is unlikely 
to be the driving source since its SED represents a reddened 
photosphere. IRS 1 1 is probably a reddened background star that 
is unrelated to the Coronet cluster. 

MHO 2007 - IRSI/HHIOOI R: MHO 2007 is a g roup of 
knots discovered in the images of lPeterson et al. I (1201 ih and our 
H2 images. Additionally, our deep H2 data reveal a faint, but 
clear bow shock feature at the same location (see Fig. 1). The 
bright knots of MHO2007 are the H? counte r parts of the HH 
objects HH 101A,B,N, and S. iPeterson et al. I 00 111) associate 
this group of knots to the red lobe originating in IRS5N. Also, 
there is a star embedded in this group of MHO features /HH 
objects that was lab e led St ar30 by Hartigan & Lada (1985). 
iHartigan & Graham I (Il987h argued that Star30 is not the 
driving source for the HH 101 group of knots. From our SED 
fitting we concur with this viewpoint: using the USNO and 
2MASS magnitudes we find that the photometry of star30 
is be st described by an extin cted photosphere rather than a 
YSO. lHartigan & Lada I ( Il985h and iHarti gan & Graham I (Il987h 
suggest that HH 101 must be part of a long coUimated outflow 
originating in HH lOOIR (close to the Coronet cluster). The 
faint, but clear, bow-shock feature detected by our H2 data could 
indeed be the terminating point of the flow from HHIOOIR. 
In Fig. 1, the tip of the arrow drawn from IRSl points to the 
center of this bow- shock feature. Thu s, the speculatio n of 
lHartigan& Graham I ([1987) and Hartiga n&"Ladal (Il985h can 
also represent a valid scenario. This flow from HH lOOIR/IRSl 
would then measure ~365", corresponding to a single lobe 
length of -0.23 pc. In such a case, HH lOOIR/IRSl may be 
driving a parsec scale outflow. 

MHO201 1 D, E - S MMlb/SMAl: The submm source 
identified by 'Nutter et al. ' (12005) as SMMlb and resolved as 
SMAl by Groppi et al. (2007) appears as a point source in 
the Spitzer-IRAC images. The luminosity (see Table. 4) of this 
source is relatively higher compared to other modeled sources in 
the region. The newly detected bow-shock features MHO 201 1, 
D & E are well aligned with this source in projection. Further, 



the PA of the proper motion vector estimated bv IPeterson et al. I 
(1201 ih points away from this source. The source SMMlc/SMA2 
appears to be aligned with MHO 201 1, but is a less probable op- 
tion because: a) we have already shown that SMA2 very likely 
drives MHO2000, b) SMAl is thought to be relatively older 
(a transitional ClassO / I object) than SMA2 (Class protostar) 
fGroppi et alTI l2007t iForbrich & Preibish I l2007h (assuming a 
similar mass), therefore expected to be associated with a longer 
and relatively fainter flow. MHO2011, D, E therefore are best 
associated with SMAl. 

MHO 2014 - SMM6: This is a bright extended H2 feature, 
with a bow-like component (to the east; lefthand side in Fig. 2) 
and a clumpy component (western half). The bow-like compo- 
nents are associated with the previously known Herbig-Haro fea- 
tures HH 732 A, B, and C (Wang et al. 2004). These HH objects 
are n ot identified with any clear driving source. IPeterson et al. I 
(I2OI 1) compute the position angle (PA) for the proper motions 
to be in the range 31°-34°+15°and associate it with the blue lobe 
originating in IRS5N. They suggest MHO2007 as the red lobe of 
this flow. This scenario may appear likely, because MHO2007, 
MHO2014 and SMM4/IRS5N align along a straight line and 
SMM4 displays a peak at millimeter wavelengths. However, the 
1 .2mm flux (and the resulting disk mass) of IRS5N, measured 
by the same authors, is weaker than sources such as SMAl 
or SMA2. Furthermore, this source is visible in the infrared 
bands and is not very luminous. Our SED modeling suggests that 
SMM4/IRS5N is a substellar object with a luminosity ~ 0.5 Lq. 
These data together do not support the scenario of IRS5N driving 
a powerful, parsec scale outflow comprising of MHO2014 and 
MHO2007. It should also be noted that P eterson et a l. I ([201 1|) 
associate SMM4 with IRS5 and state that the new peak detected 
using SMA observations is associated with IRS5N, while SMM4 
and the new SMA peak are indeed the same, coinciding with 
IRS5N. 

Here we propose SMM6 (MMSlO-11) as an alternative 
drivin g source for MHO2 014 dNutteret all 120051: IChini et al. I 
'2003|). IChini et al. I d2003l) note that there are n o near-IR or 
FIR s ources associated with M MSlO-11 peaks. dNutter et al. I 
120051) discard the proposition of fChini et al. I d2003l) that MMS6 
is a deeply embedded protostar based on their interpretation 
of this source as a prestellar object. This interpretation may 
be valid, since the Spitzer MIPS 24jum or 70jum images do not 
reveal point-like or extended source associated with SMM6. 
However, the 160/im emission in this region is bright and satu- 
rated, rou^hlyToUowingthe morphology of the 850yum emission 
from iNutter et all d2005l) . Our supposition that SMM6 is the 
driving source for MHO2014 requires a PA of 4 1°, which is wel l 
within the value and error bars quoted by Pe terson et al. I d201 lb . 
Given that this is one of the brightest MHO in the region con- 
forms to the proposition that it should be driven by a deeply em- 
bedded protostellar object invisible at the infrared wavelengths. 

The western half (righthand side in Fig. 2) of MHO 2014 is 
very clumpy and not likely to be driven by a well-collimated 
ClassO jet/outflow. The morphology of this clump and its 
location in the vicinity of the bright young stars TY CrA 
and HD 176386 may imply a different origin. TY CrA is a 
well known variable star of Orion type, and HD176386 is a 
pre-main-sequence double star coinciding with a submm peak 
dChini et al. II2003I) . Lacking proper motion data, we therefore 
speculate that this clumpy knot might be due to an unidentified 
component of the TY CrA or due to the double components of 
HD176386. 
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30.0 19:02:00.0 30.0 01:00.0 

RA (J2000) 

Fig. 2. Continuum-subtracted H2 image of the mapped region, displayed with Hnear gray scales. Five groups of H2 line emission 
features are marked by ellipses. 



Table 2. MHO and source identifications and flow lengths 



MHO 


RA (J2000) 


Dec (J2000) 


Length 


Source 


Comments 


2000 A, B, C 


19 02 04.8 


-36 54 45 


198" 


SMA2 


HH99 A, B, C 


2001 


19 02 01.00 


-36 56 38 


113"/145" 


IRS6 or IRS5N 




2002 


19 01 59.1 


-36 57 17 


70" 


RCrA or CX034 




2003,2007,2008 


19 01 46.0 


-37 00 06 


350" 


IRSI/HHIOOIR 


HH96, HH97, HHIOO, HHlOl 


2004 


19 0146.8 


-36 56 55 


64" 


IRS6 and/or RCrA 


HH104(jet like) 


2005 A, B 


19 01 25.9 


-36 57 20 


258" 


IRS5a/b 


thin faint arc 


2005,C,D,E.. 


19 01 15.6 


-36 57 25 


87" 


SCrA 


HH82 A, B 


2006 


19 01 39.2 


-36 59 02 


330" 


IRS2 


HH730 A, B, C 


2008 D 


19 01 40.0 


-37 00 05 


16" 


Star 1-100 


1-lOOHH 


2009 


19 02 12.2 


-36 56 43 


193" 


IRS7A 


HH735, HH734 


2011 D,E 


19 02 12.2 


-36 56 43 


193" 


SMA1/IRS7B 


HH735, HH734 


2011 A,B, 2012 


19 02 20.6 


-36 55 55 


275" 


SMM2 


faint bow shocks 


2013 


19 02 09.2 


-36 54 53.0 


289" 


TCrA 


northern lobe, HH733 (bowshock) 


2015 


19 01 55.7 


-36 58 42 


67" 


TCrA 


southern lobe bowshock 


2014 


19 01 57.5 


-36 52 12 


236" 


SMM6A1MS10-11 ? 


HH732 


2016 


19 01 39.5 


-36 53 02 




TYCrAandHD176356 





5. Discussion 

We discovered new MHO features tracing individual flows in 
the northern part of the Corona Australis region, which encom- 
passes the Coronet cluster Together with the previously known 
flow components, we have cataloged 16 individual flow lobes 
in the studied region (Table. 1). Eleven of these flows are local- 
ized in the Coronet cluster region, which has a membership of 
about 14 stars that are detected in the x-ray dForbrich & Preibishl 



12007 ^ near-IR, mid-IR, FIR (Spitzer) ([P eterson e t all 1201 lb 
and submm photometric da ta dNutter et al. i,2005; .Groppi et al. I 
I2OO7I: iPeterson et £1711201 ih . Thus, roughly 80% or more of the 
embedded sources in the Coronet cluster display an active out- 
flow component. The region m apped in H2 here is ro u ghly s im- 
ilar to the region covered by iForbrich & Preibish ] (I2007h m 
the x-ray and Spitzer data. These authors identify nine embed- 
ded objects (meaning Class 0/1 and flat-spectrum sources) in 




Fig. 3. Continuum-subtracted H2 image of the Coronet subregion dis played using a logar i thmic scaling. The circular syrnbols sh ow all sources detected by the Spitzer IRAC 
images. The identifications are reproduced from WiMngeHU dHHI) (IRS), iGroppi et all (|2007|) (SMA), and lNutter et al. I (|2005|) (SMM). The aiTows join the pr obable driving 
source with the corresponding MHO. When two arrows originate in the same source, it implies both outflow lobes are visible. The CO bipolar outflow mapped bv lGroppi et al. I 
(12004) is shown by solid (blue lobe) and dotted (red lobe) contours. For the sake of clarity, MHO 2000, which is situated close to SMM2, is not marked in this figure. 
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Fig. 4. The observed SEDs and their best-fit models. Top left panel shows an example of an extincted photosphere while the rest are 
YSOs. Circular symbols represent the data points, upward-facing triangles are lower limits and downward facing triangles are upper 
limits. The dashed line corresponds to the stellar photosphere model used in the best YSO mod el fit. The gray li nes show the models 
that satisfy the;\f^ criteria described in the text. For the source IRS5N, the SMA data bv lPeterson et al. I (1201 Ih is additionally used 
in the SED fitting. 



this region. Together with the newly ide ntified SMA sources 
(iGroppi et al. Il2007t iPeterson et 31711201 ih . all of the embedded 
sources are associated with MHO flows. 

If we compare the submm cores from iNutter et al. I (l2005h 
with the detected MHO flows in this region, all cores except 
SMMIA and SMM6 ai-e associated with a detected infrared 
source driving an outflow. SMMIA and SMM6 are the two cores 
that these authors classified as prestellar, of which we associate 



SMM6 with MHO2014, even though no infrared source is de- 
tected with this core up to 24jum . Therefore, the studied region 
contains one, or at most two, starless cores and eight protostellar 
cores, which implies roughly 10% of starless cores in the re- 
gion. The Corona Australis core is therefore a good example of 
a young star-forming cluster, with a very high rate of detected 
outflow activity and young stars. 
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Table 3. SED-fitted parameters of selected sources in the Coronet. For each parameter, the minimum, maximum, and weighted 
mean (inside parantheses), selected by the;i'^ criteria discussed in the text are displayed. These values represent the group of models 
shown with gray lines in Fig. 4. SMM ident i ficatio ns are from Nutter et al. (2005), I RS identifications fr om Wil king et al. , (.19971) . 
SMA identifications are from lGroppi et al. I (l2007l) . and CXO identification is from lForbrich & Preibishl (l2007h . 



source 


Class 


L 


M. 




Mrf,-,i. 






log(Lo) 


Mo 


log(Mo>'f--') 


log(Mo ) 


IRS5N/SMM4 


o/i 


-0.2, 1.8, (0.2) 


0.1,2.9, (0.5) 


-8.9, -6.0, (-7.3) 


-3.8,-1.6, (-2.2) 


IRS7A/7W 


o/i 


0.1,2.8, (1.1) 


0.1, 6.8, (1.2) 


-9.9, -4.7, (-7.0) 


-4.1, -0.4, (-2.4) 


CX034 


O/I 


-0.4, 2.5, (0.4) 


0.1,6.5, (0.5) 


-11.8, -5.2, (-8.0) 


-5.0, -0.8, (-2.8) 


SMM 


O/I 


1.0, 2.9,(1.9) 


0.2, 8.0, (3.3) 


-9, -4.6, (-6.9) 


-3.8, -0.5, (-1.7) 


SMM2 


O/I 


-0.3, 2.0, (0.4) 


0.1,4.2, (0.7) 


-10.5,-6.4, (-8.1) 


-4.4, -1.3, (-2.8) 



From Fig. 3, we can see that for each driving source, only 
one of the outflow lobes is traced either through H2 emission 
or HH objects. T CrA is the only exception in this region, as- 
sociated with two lobes identified with bow shock features. For 
all the other flows, the identified flow lobes are oriented in such 
a way that they are flowing in a direction away from the cen- 
tral dense core and entering the less dense medium. The out- 
flow lobes flowing towards the central dense core is usually 
absent. The dense core in the Coronet cluster is rough ly cen- 
tered on the Herbig Ae/Be star RCrA (IChini et al. ''2003). We 
can see from Fig. 3 that the flow lobes from the sources to the 
left (east) of RCrA are all moving towards the east or northeast, 
while their eastern/southeastern counter-lobes are not detected. 
The sources located to the right (west) of RCrA display a simi- 
lar, complimentary behavior. It appears that the outflows driven 
by the young stars are "stunted" because of the increased power 
needed to penetrate the dense gas inside the core. In compari- 
son, the flow from TCrA, which is orientated north-south and 
avoids the dense core, is bipolar The dense core around Coronet 
is concentrated within a 3'-4'diameter. Since the detected single 
lobe lengths surpass this value, even assuming different inclina- 
tions, the brightest bow shocks or knots from the counterlobes 
should have been detected in our deep H2 images, if the "stunt- 
ing" effect was not real. Comparing Coronet to NGC1333, we 
can see further support of the likely "stunting" effect scenario. 
NGC1333 is anot her nearby region w ith a high incidence of de- 
tected MHOs (see lDavis et al. 1120081 and references therein). In 
NGC1333, MHOs are associated with bipolar lobes of multi- 
ple ou tflows, even in a relatively shallow survey of Davis et al. 
( I2OO8I) . However, the dense cores are ordered along a curved 
filament with gap s of low density regions (see Figs. 3 & 7 of 
IWalsh et al. ll2007h . and the bipolar outflows can be seen punch- 
ing out roughly perpendicular to the filamentary axis. The MHOs 
detected in this region can t hen be verified by comparing with 
Fig.5 of iDavis et al. ] dloosh . Therefore, we speculate that the 
higher extinction dense core, has an important effect on the opti- 
cal and infrared detection of the flow lobes. Only the lobes mov- 
ing away to the rarer medium produce significant H2 emission 
resulting in good detections. This observational feature may be 
a qualitative indication of the lower density of the outflowing gas 
compared to the density of the natal dense cores in which they 
form. It may also demonstrate how the outflowing gas penetrates 
out into a region of least pressure. 

The nature of driving sources and associated flows: It can 
be seen from Table 2 that more than half of the detected flows, 
in particular the flows originating from the Coronet, display a 
single lobe length of >200". This implies projected total flow 
lengths of >0.1pc for the bipolar outflows at a distance of 130pc. 
By assuming roughly similar inclination angles to the observed 
flows, a canonical jet speed of 100 km s ', a flow length of O.lpc 



indicates a dynamical time scale of ~1200yr for the outflows. 
Since the molecular hydrogen flows are driven by YSOs in their 
early evolutionary stage and since roughly 80% of such identi- 
fied members are located in the Coronet cluster and are driving 
H2 flows, the cluster should be very young with an age -IMyr 
Mever & Wilking (2009) have evaluated the masses and ages of 
the Coronet sources using near-infrared spectra and evolutionary 
models. They find an age spread of 0.3-3 Myr and masses in the 
range O.2-2.5M0for the young aggregate. 
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Fig. 1. A coloured image of the region studied in this work. The colour image was obtained by coding Spitzer IRAC 4.5yum image as red, H2 narrowband image as green and SII 
image (IWang et al. II2004I) as blue. The SII data do not cover the lower portion of the image. 



